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Non-Gaussian morphology on large scales: 
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Abstract. In order to quantify higher-order correlations of the galaxy cluster distribution we use a complete 
family of additive measures which give scale-dependent morphological information. Minkowski functionals can 
be expressed analytically in terms of integrals of n-point correlation functions. They can be compared with 
measured Minkowski functionals of volume limited samples extracted from the Reflex survey. We find significant 
non-Gaussian features in the large-scale spatial distribution of galaxy clusters. A Gauss-Poisson process can be 
excluded as a viable model for the distribution of galaxy clusters at the significance level of 95%. 
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1. Introduction 



Galaxies: clusters: general - Cosmology: observation - Cosmology: 



tant leading to a non-Gaussian density field. Also the 
process of galaxy formation may introduce n on-Gaussian 



The s patial distribution of galaxy clusters poses impor- featu res if the "biasing" is non-linear (see e.g. |Scoccima7ro 



tant cdnstramts on cosmologicai models. The abundance 
of clusters and especially its evolution with redshitt is very 
sens itive to parameters of t he cosmologicai models (see 
e.g. Kitayama fc Suto 1997, Borgani et al. 1999, Bahcall 



2000 ). Typically one argues that on large scales, the evo- 



200C, Kerscher et al. 2001a). To quantify the large-scale 



structures traced by the galaxy clusters we have to go be- 
yond the number density. 

Scenarios describing the formation of structures in the 
Universe start with a mass density field showing only small 
deviations from the mean density. Inflationary scenarios 
suggest that these density fluctuations can be modeled 
as a Gaussian random field completely specified by its 
mean value and the power spectrum or two-point cor- 



lution is still in the linear regime, and one expects that 
the smoothed density field is proportional to the initial 
Gaussian field. However, structures like walls and fila- 
ments were observed in the galaxy distribution on large 
scales ( Huchra et al., 199C| ; Shectman et al., 1996 ). These 
non-Gaussian features appear at a low density contrast 
and are therefore hard to detect. The sensitivity of the 
Minkowski functionals, even if only a small number of 
points is available, allows us to quantify the non-Gaussian 
morphology of these structure on large scales. Walls and 
filaments were predicted by analytical and numerical work 



based on the Zel'dovich approximation (Zcl'dovich, 1970 



relation function (e.g. Kolb & Turner 1990). In the ini- Arnol'd et al., 1982; Doroshkcvich et al., 1996) and re- 



tial stages of structure formation the linear approxima- lated approximations (Kofman et al., 1992; Bond et al. 
tion is often used to evolve these fiuctuations preserv- 1996). A^-body simulations could verify that these struc- 



ing their Gau ssian nature a nd increasing their amplitude 
only (see e.g. Peebles 1980 ). With growing over-density 
the nonlinear couplings become more and more impor- 



tures are generic features of the gravitational collapse for 
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Cold Dark Matter ( CDM) like initial c onditions ( [Melott 
& Shai |idarin, 1990| ; [Jenkins ct al., 19981 ). 

Since observations supply us with the positions of 
galaxies and galaxy clusters in space, our methods will 
use this point distribution directly. No smoothing is in- 
volved. Therefore we have to give a clear definition what 
a "Gaussian" point distribution, the Gauss-Poisson pro- 
cess, is. Some of the statistical properties of random fields 
directly translate to similar statistical properties of point 
distributions, but also important differences show up. The 
equivalence of the Gauss-Poisson process with a simple 
Poisson cluster proc ess, allows us to simulate a 'Gaussian" 
point distribution ( Kerscher, 2001 ). With these simula- 
tions we will perform a Monte-Carlo test to determine 
the significance of the non-Gaussian features in the clus- 
ter distribution. 



discipline of integral geometry (see e.g. Hadwiger 1957) 
supplies us with a suitable family of such descriptors, 
known as Minkowski functionals. These measures embody 
information from every order of the correlation functions, 
are numerically robust even for small samples, and yield 
global as well as local morphological information. The 
Minkowski functionals are additive measures which allows 
us to calculate them efficiently by summing up their lo- 
cal contributions although they depend on all orders of 
correlation functions. The application of Minkowski func- 
tio nals in statis tical physics and cosm ology are reviewed 
by iMecke (200"o| ) and [Kerscher (2000D , respectively. 

Samples of galaxy clusters are based mainly on opti- 
cal observations, where the clusters are selected as galaxy 
over-densities in the two-dimensional maps on the celes- 



Sta mistical measures provide important tools for the 
comparison of the large-scale structure in the Universe 
with theoretical models. The discriminative power of this 



tial s pher e (c.f. [Abell 1958| , \Ahell et al. 1989[ [Palton et al 
1997, and Gal et al. 200C). Projection effects seem to have 



a non-negligible effect on the sta tistical analysis of these 
optically selected cluster samples ( Katgert et al., 1996 ; van 



compa lson depends chiefly on the statistical measure. Haarlem et al., 1997). Only in recent years X-ray selected 



The most frequently employed measure was and still is 
the two-point correlation function, or the power spec- 
trum. Both are nowadays an imperative in the analysis 
of any galaxy or cluster catalogue: for the Reflex clus- 
ter catalogue see [Collins et al. (2000 ) and Schuecker et al 
(200l|2 They give important information on the fluctu- 
ation spectrum of matter. However, they appear to be 
blind to morphological features. Indeed, completely differ- 
ent spatial patterns and point distributions could display 
the same two-point correlation function, i.e., no direct 
conclusions about the morphology of the structure can be 
drawn from an analysis with these two-point measures 



( [Baddeley fc Silverman, 1984[ ; [Szalay, 1997[ ; [Pan fc Coles ^ 
200 q; [$:erscher, 2001[ ). Highe r-order correlation functions 
immediately come to mind if one wants to go beyond the 
two-point correlation function. And indeed three-point 
correlati ons were detected in the distribution of galaxy 
clusters ([Toth et al, 1989[) . However, there is a concep- 
tual problem since n-point functions {n > 3) depend on 
3(n— 1) — 3 parameters even for isotropic and homogeneous 
point distributions. Already for the three-point correla- 
tion function we are not aware of a study where its depen- 
dence on all three parameters was estimated. Clearly, in- 
tegral information is mandatory and necessary. This may 
be accomplished e.g. for the three-point function by av- 
eraging over the shape of triangles, or by considering the 
(factorial) moments of counts in cells (see e.g. Peebles 1980 
and Szapudi & Szalay 1993). Another effort to go beyond 
the two-point correlation function comprises the perco- 



cluster samples have been completed. Since the X-ray lu- 
minosity is proportional to the baryonic density squared, 
over-densities are more emphasized. Consequently, the 
contamination of the catalogue by chance alignments 
due to projections is reduced (Bohringer et al., 2001). 
Assuming a virial relation, the X-ray luminosity of the 
galaxy cluster can be related to its mass. 

2. Morphology of large scale structure 

Minkowski functionals have been introduced to cosmol- 
ogy as a tool to quantify the m orphology of large-scale 
structures by Mecke et al. (1994 ) where also a flrst anal- 
ysis of the distribution of galaxy clusters based on the 



Abell ct al. (198£ ) sa mple with a redshift compilation by 



Postman et al. (1992 ) was presented 



With Minkowski functionals we quantify the morphol- 
ogy of a sufficiently well behaved compact body K C M.'^ 
by assigning it a number M^{K) E R. The Minkowski 
functionals are motion invariant 



M,{gK) = AUK), 



(1) 



lation analysis (3handarin, 1983). Also the genus, closely 
related to the Euler characteristic, is often employed to 
quantify deviations from a Gaussian density fields (see e.g. 



Hamilton et al. 1986, Melott 1990 and references therein) 



where 5 = (x, 0) are the movements in three dimensions, 
i.e. translations x and rotations Q. As already emphasized 
the additivity property 

M^{K U K') = AU{K) + A'U{K) - A'U{K n K') (2) 

serves as the construction principle of these measures. All 
the Minkowski functionals have a straightforward inter- 
pretation in terms of geometrical and topological quanti- 
ties as summarized in Table |l|. Minkowski functionals are 
distinguished from other geometric measures by the theo- 



For the construction of statistical methods, sensitive to 
the large-scale structures, additivity is a heuristic princi- 
ple which can guide us how to define useful measures which 
do not depend on all these parameters. Additivity yields 
robust, local decomposable measures. The mathematical 



rem of Hadwiger (1957), which states that there are only 
four independent scalar functionals in three-dimensional 
space, which are motion invariant, additive, and continu- 
ous for convex bodies. Hence, every additive geometrical 
measure M{K) which does not depend on the position and 
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Table 1. Minkowski functionals and the normalized 
Minkowski functionals <f>^ (Eq. (||)) in three-dimensional 
space expressed in terms of more familiar geometric quan- 
tities. 



geometric quantity 


M 






volume 


V 





V 


l//(^r^iV) 


surface area 


A 


1 


A/8 




integral mean curvature 


H 


2 




H/{A-KrN) 


Euler characteristic 


X 


3 


3x/(47r) 


X/N 




Universe. Since there are only four numbers of the MFs 
in three dimensions, compared to a correlation function, 
it is necessary to define morphological functions M^{r). 
The proper technique to do this for general random spa- 
tial structures are (erosion) dilation operations (see Fig. ^ 
and 3erra 1982). In case of point patterns this tech- 
niques reduces to fixing balls Br of radius r at each point. 
With Minkowski functionals we quantify the geometry and 
topology of union set of these balls Ar = UiLi ^'•(xj). 
The radius r is employed as a diagnostic parameter. In 
such a way, we obtain scale-dependent integral infor- 
mation on higher-order correlations of the distribution 
of galaxy clusters and not only two-point correlation 
functions of the large-scale distribution. Erosion/dilation 
techniques combined with additive Minkowski functionals 
have been successfully applied in many areas, including 
condensed m atter physics ( Mecke, 2000 ), g eology (Arns 
et al., 200l4 |b|), and digital image analysis ( Serra, 1982 , 



Fig. 1. Dilated points: spheres of varying radius attached 
to the galaxy cluster from the Reflex sample L12 (see 
Table H). 

orientation of the body K in three-dimensional space, can 
be written as a linear combination of the four Minkowski 
functionals: 



Af(i^) = ^c,Af,(i^). 



(3) 



i/=0 



Beisbart et al. (2001| ) discuss the vector valued extensions 
of the Minkowski functionals. 

The cluster distribution provides us with a point set 
X = {xijf^j in three-dimensional space. One may think 
of X as a skeleton of the large-scale structures in the 



2.1. The Boolean grain model 

The simplest model of a random point distribution is the 
Poisson process. By attaching grains to each of the points, 
in our case balls of radius r, we arrive at the Boolean 



grain model. Mecke fc Wagner (1991 ) presented a method 
to calculate mean volume densities of the Minkowski func- 
tionals for this model. We repeat their arguments, since its 
extension allows us to calculate the Minkowski functionals 



of correlated grains in Sect. 3.1 



Iterating the additivity relation (Q) one obtains for the 
union Ar = [JiLi Bri'x.i) of N spheres = Br{'X-i) of 
radius r and center x, the inclusion-exclusion formula 



M,{Ar) = ^^-(B^) - H ^MB^ n Bj) + 

i i<j 

+ ...i-if+'AUiB,n...nBN). (4) 



Generally, a point process in a domain Q with vol- 
ume |r2| is specified by a sequence of product densities 
p„(xi, ... , x„) with the mean number density g = gi{xi). 
Qni'X-i, . . . ,x„)dVi . . . dVn is the probability of finding n 
points in the volume elements dVi to dVn. The volume 
density m,^{Ar) of the v-th Minkowski functional per unit 
volume for the augmented coverage Ar are then obtained 
from the inclusion-exclusion formula (^) in the form 

m^(A; {gn}) = 

J "^^^ [ n j SniTn) (5) 



where we introduced, for convenience, the variable r„ = 
(xi,... ,x„) with the integration measure J^dTn = 
YVi^i Jfi'^^i- Obviously, the Minkowski functionals em- 
body information from every order n of the n- 
point densities p„(xi,... ,x„). If the product densi- 
ties gn(xi,... ,x„) — were independent of position 
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(Poisson distribution of density g = N the integrals 
in Eq. (j^) can be performed using the fundamental kine- 
matic formula ( Blaschke 193(^ , Santalo 1976| ) 



100 



LOS L12 L20 L30 



J dg M,{K n gK') - ^ Q M^{K)M,^^{K'). (6) ^ 



Equation (||) describes the factorization of the Minkowski 
functionals of the intersection KOK' of two bodies K and 
K' if one integrates over the motions g = (x, Q), i.e. trans- 
lations X and rotations Q oi K'. For Poisson distributed 
spheres Br of radius r one obtains the mean values of the 



Minkowski functionals per unit volume (Mecke & Wagner 
19911, Btoyan et al. 1995[) 



in^{Ar]g) = ^< 1 - exp 



t=0 

(7) 



and particularly in three dimensions the normalized 
Minkowski functionals 



miy{Ar) 
M,{Br) Q 



(8) 



read with rj = g 



(9) 



Mr\Q) = (l-e-")/?7 , 
^i{r;g) = c"'' , 
$2(r;f?) = (l-2|i^)e-') , 
^3{r;g) = [1 - 3v + ^V^^-" ■ 

2.2. The Reflex duster sample 

The construction of the Reflex cluster sample is de- 
scribed in detail by Bohringer et al. (2001 ) and the statis- 
tics of the cluster distribution is described by Collins 
et al. (200C ) for the tw o-point correlation function and in 
Schuecker et al. (2001) for the density fluctuation power 
spectrum. The survey area covers the southern sky up to 
the declination 6 < 2.5° avoiding the band of the Milky 
Way, \bii\ < 20° and the regions of the Magellanic clouds. 
The total survey area is 13924 deg^ or 4.24 ster. Tests in- 
cluding the number counts (log iV log S'-function) , the co- 
moving densities, (V/Vmax) tests, and comparisons to sim- 
ulations, described in the above mentioned papers, show 
that the selection function is well documented. 

The X-ray detection of the clusters is based on the 
second processing o f the RASS (ROSAT All Sky Survey, 
Voges et al. 1999) exploiting a primary (MPE inter- 
nal) source detection list comprising 54076 sources in the 
Reflex area down to a detection likelihood of L > 7 (see 



Voges et al. 1999). For all these sources the X-ray param- 



eters are reanaly sed by the growth curve analysis method 
as described by Bohringer et al. (2000| ) which provides 
a flux measurement with significantly less discrimination 
against extended X-ray sources than provided by the stan- 
dard analysis of the RASS. The results of this reanalysis 




0.001 



Fig. 2. Luminosity - redshift distribution of Reflex 
clusters of galaxies (points) and the applied ranges for 
the extraction of volume-limited subsamples. The redshift 
and luminosity intervals of the respective volume-limited 
subsamples L05, L12, L20, and L30 are marked by con- 
tinuous, short-dashed, long-dashed, and dotted lines. The 
subsamples are described in Table ^. Note that for con- 
ventional reasons the X-ray luminosities are given in units 
of iJo = 50kms-iMpc"\ 



are used to produce a flux-limited sample of RASS sources 
with a nominal flux > 3 • 10^^^ erg s^^ cm^^ in the 
energy band (0.1-2.4 keV). 

The cluster candidates are finally identified or removed 
from the sample as non-cluster sources by a detailed 
documentation of the X-ray and optical source proper- 
ties, literature information, and spectroscopic informa- 
tion including redshift measurements obtained by follow- 
up observations within the frame of an ESO key pro- 
gram. Further tests of the sample completeness based on 
a search for clusters among the significantly extended X- 
ray sources and a search for X-ray emission from clusters 



cataloged by Abell et al. (1989) independent of the RASS 
source detection supports the completeness estimate of 
> 90% for a flux-hmit of 3 • 10~^^ergs"^cm~2. The high 
completeness concerning the optical identification makes 
the data set an effectively X-ray selected sample of galaxy 
clusters. The final cluster sample includes 452 clusters and 
there are three objects left in the list with uncertain iden- 
tifications and redshifts. These three objects are excluded 
here in the further analysis. 

For the determination of MFs complete volume- 
limited subsamples are needed. The Reflex cluster sam- 
ple is per construction X-ray flux-limited so that the 
fraction of luminous clusters increases with redshift (see 
Fig. 1^). Volume-limited distributions are selected by intro- 
ducing upper redshift and lower luminosity limits (vertical 
and horizontal lines in Fig. |^). In order to reduce possible 
(error mig ration) effects wh ich might occur at the flux- 
limit (e.g., Eddington 1940| ) the upper redshift limits are 
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Table 2. The volume-limited cluster samples used in our 
investigations consisting out of N clusters closer than R 
with X-ray luminosity higher than Lmin- The constraint 
gC2 will be explained in Sect. 0. 



Sample 


R 




N 




gC2 




[Mpc/h] 


[erg/s] 








L05 


180 


0.5 X 10** 


74 


8.9 X 10"'' 


3.2 


L12 


260 


1.2 X lO^-* 


95 


3.8 X lO"'* 


2.2 


L20 


330 


2.0 X 10'''* 


86 


1.7 X 10"'* 


0.6 


L30 


385 


3.0 X lO""* 


62 


0.8 X 10"^ 


0.3 



set slightly below the formal redshift limit, especially for 
large luminosities where the effects could be largest. 

The completeness of the differ ent volume-limited sub- 
samples is illustrated in Fig. 3 in Bchuecker et al. (2001 ) 
which includes the subsamples denoted by L05 to L30 
in Table 0. Similar volume-li mited samples as lis ted 
in Table g have been used by Collins et al. (200C ) in 



their analysis using the two-point correlation function. 
Comoving distances have been calculated according to the 
Mattig formula with = 1, /i = 0.5, and A = 0. The 
flat redshift-independent distribution of comoving clus- 
ter number density suggests the absence of large incom- 
pleteness effects of the subsamples in the redshift range 
< z < 0.15. We thus expect no significant artificial fluc- 
tuations introduced by incompleteness effects on scales up 
to comoving radial distances of i? < 400 /i"^ Mpc. 

2.3. Minkowski functionals of tlie Reflex clusters 

To study the morphology of the large-scale distribution of 
galaxy clusters we consider a series of volume-limited sam- 



ples from the Reflex cluster catalogue ( Bohringer et al 
2001). The volume densities 



of the Minkowski func- 



tionals were calculated using the minus-sampling bound- 
ary correction, based on partial Minkowski functionals 
as suggested by [Mecke 
Schmalzing 



et al. (1994) (for details see 
Diaferio 200C| ). The survey is bounded by 



6 < 2.5° and |6| > 20°, but also several regions in the 
small and large Magelanic clouds were excluded from the 



sample (for details see Bohringer et al. 2001). To estimate 
their influence on the Minkowski functionals of the sam- 
ples we filled these regions with random points with the 
same number density. The comparison in Figs. shows 
nearly identical results for filled or unfilled regions in the 
Magelanic clouds. 

The overall features seen in the Minkowski functionals 
of the Reflex clusters are similar to the one observed 



in the Abell/ACO cluster sample of Plionis & Valdarnini 
|(199l| ) as analyzed by Kerscher et al. (1997). We are lim- 



ited by the smaller sample size and the boundary cor- 
rection used. Only few galaxy clusters contribute to the 
Minkowski functionals for large radii. Therefore we are 
not able to trace the large-scale structure to the limit 
where the sample volume is filled by the union set of 



0.3 



0.2 





5 10 15 20 
r [h iMpc] 



5 10 15 20 
r [h 'Mpc] 



2.5xl0-s 




5 10 15 20 
r [h-'Mpc] 



5 10 15 20 
r [h"'Mpc| 



Fig. 3. Minkowski functionals of the volume-limited sam- 
ple L05 (solid line) compared to the Minkowski function- 
als of a Poisson process with the same number density 
(dotted line, gray shaded one-cr area). The numerically 
determined mean is in perfect accordance with Eq. (^. 
The dashed one-tr area is obtained by filling the excluded 
area around the small and large Magelanic clouds with 
randomly distributed points of the same number density 
as the rest of the sample. 




Fig. 4. Minkowski functionals of the 
pie L12. Same conventions as in Fig. 



volume-limited sam- 
I 



6 



M. Kerscher et al.: Non-Gaussian morphology from the Reflex cluster catalogue 




20 30 
[h-iMpc] 




10 20 30 
r [h^'Mpc] 



40 



Fig. 5. Minkowski functionals of the volume-limited sam- 
ple L20. Same conventions as in Fig. ^. 
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a clear trend from strong clustering towards only small 
differences from the Poisson distribution. Increasing the 
depth of the volume-limited samples we allow for clus- 
ters with higher X-ray luminosity. Considering the ampli- 
tude of the two-point correlation function, galaxy clusters 
with higher X-ray luminosity should show stronger cor- 
relations (Kaiser, 1984; Bardeen et al., 1986). However, 
in the deeper volume-limited samples, with the more lu- 
minous clusters, also the number density decreases. The 
sparseness of the point distribution competes with the in- 
creased amplitude of the two-point correlation function. 
Indeed, for quite general conditions, a point distribution 
converges towards a Poisson process under thinning, i.e. 
unde r randomly deleting points (e.g. Daley & Vere-Jones 



1988, Sect. 9.3). For the Minkowski functionals the behav- 



ior in sparse samples may be explained by considering the 
expansion of the normalized Minkowski function als <&„ in 
terms of 77 = g(47r/3)r^ around zero (see also Kerscher 
et al. 2001b ). Based on the expansion (||) of MFs in terms 



Fig. 6. Minkowski functionals of the volume-limited sam- 
ple L30. Same conventions as in Fig. ^. 



balls. Additionally to the Minkowski functionals of the 
clusters the results for a Poisson process with the same 
number density inside the sample geometry is shown in 
the Figs. H-|6[ Increasing the depth of the volume-limited 
samples from L05 to L30 the Minkowski functionals show 



of n-point densities one gets to the lowest order in r] 



*m(A) = 



2r 



ds 



,{r;s)[l + Us))+0{7j'). (10) 



The functions 



h{r;s)= 1-^, 



lf^)3 
2\2r) ' 



/2(r;s)- 1 
/3(r;s) = e(2r-s) 



2r 



(11) 



resemble the Minkowski functionals I^ir; |x|) = 
Mf,{BriO) n Brix))/Mf,{Br) of the intersection of 
two spheres of radius r and distance |x|. 0{q) is the 
step function equal zero for q < and one for q > 0. In 
Eq. ( p^ the integral and its pre-factor give a dimen- 
sionlcss geometric number depending on the correlation 
function ^2(5) and the measure Mfj_. Terms proportional 
to 77" include intersections of n spheres weighted by the 
n-point densities. Only the two-point correlation function 
is important for small 77, but the Minkowski functionals 
become increasingly more sensitive to higher-order 
correlations with larger rj. For very small r/ we essentially 
arrive at a Poisson process as numerically verified by 
Kerscher et al. (1998 ). A small rj may be obtained either 



with a small radius r of the spheres, or, as in our case, by 
a low number density g. 

Now let us describe the features in the MFs in more 
detail. The strong clumping in the distribution of galaxy 
clusters is causing the lowered values of the volume den- 
sities rrif^ of the Minkowski functionals compared to the 
Poisson values. In a clustered point distribution, the 
spheres in the union set Ar = Ui!=i Br{^i) overlap signifi- 
cantly already for small radii. This is leading to a reduced 
density of the volume ttiq, surface area ttt-i, and integral 
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mean curvature 7712. The density of the Euler characteris- 
tic ma decreases since for small radii mainly the number 
of connected objects is counted - no tunnels and cavities 
have formed yet. A tunnel through the body Ar gives a 
negative contribution of minus one to the Euler charac- 
teristic. In the sample L20 we observe the zero crossing of 
the Euler characteristic indicating that an interconnected 
network of tunnels, a sponge-like, bi-continuous "cosmic 
web" has formed for radii around 35ft-^^Mpc. 

In the deeper samples L20 and L30 the volume density 
niQ shows a tendency towards increased values compared 
to a Poisson process. With our estimator for the MFs, we 
successively shrink the sample proportional to r, where 
r is the radius of the spheres Br- Therefore, we mainly 
probe the ceuLi'al I'cgioii of Lhc sample foi' lai'gc r. The in- 
creased mo iH cauHcd by gi'adieiitH lii the numbei' dciiHlty q 
of the Reflex cluster sample, specifically the local under- 
density of clusters out to approximately lOO/i^^Mpc (see 



with r = |x2 — xi|, and in general 



There is no easy relation between the scale s of fluctu- 
ations in the number density as probed by £,2(3) and the 
radius r of the spheres used in the MF analysis. As can 
be seen from Eqs. (||) and (|l^) weighted integrals over all 
scales contribute to the MFs at a given radius r. However, 
with the radius r of the spheres Br we probe the geom- 
etry and topology of the cluster distribution in a scale- 
dependent way. The radius r can be regarded as a geomet- 
rical scale, e.g. the radius of the first zero of the Euler 
characteristic m^^Arj,) = is an estimate of the perco- 
lation threshold for our system of mono-disperse spheres 
( Mecke fc Wagner, 1991| ). At this scale Vp the large-scale 
structure elements (walls, filaments, clusters) form a per- 
colating network. 



3. Higher-order correlations in point distribution 

3.1. Minkowski functionals for correlated grains 



As for the Poisson process (Mecke & Wagner 1991, and 
Sect. 2.1) one may calculate the Minkowski functionals 
of correlated grains, in our case spheres centered on a 
clustered point set (Mecke 1994, Schmalzing et al. 1999| , 
Mecke 2000| ). Our main analytic result Eqs. (^J) and (|16[) 
expresses the Minkowski functions (r) in terms of cen- 
tered correlation functions which allows a direct compar- 
ison of measured functions with a Gaussian model where 
higher correlations are set to zero. 

The expression ^ may be used to calculate the 
Minkowski functionals for correlated grains given the n- 
point densities gi„(xi,... ,x„) of the point distribution. 
An alternative and sometimes more convenient expression 
for the densities m^{Ar) than Eq. (||) can be obtained in 
terms of the cumulants, the connected or centered correla- 
tion functions £,n{^n) with ^i(xi) = 1. For the two-point 
correlation function we have 

6(0 + l=6(xi,X2) + l = ^%^, (12) 



\-p\ 

{V} 1=1 



(13) 



Thus, the product densities Qn of order n is given by a 
sum over all possible partitions V of the coordinates r„ — 
(xi,... ,x„) into |P| parts of pi elements. Each vector 
Xj e r„ occurs exactly once as an argument x^- S Tp. of a 

cumulant Cpii^Pi) on the right side, i.e., X]i=iPi = 

Using the additivity relation (^) and the kinematic for- 
mula (||) of the Minkowski functionals one can follow the 
derivation in Mecke fc Wagner (199l[ ) so that one imme 



diately obtains the expression for the intensities (Mecke 
1994|) 



Schuecker et al. 2001, and for galaxies Zucca et al. 1997). _^ 



1 — exp 



t=0 

(14) 

due to the factorization of the integral in Eq. (^ (compare 
with Eq. (0)). The normalized specific Minkowski func- 
tionals 

Af.(r;e,U4) = E^^^^ 
^ — ' n! 

n=l 

n 

describe the deviation from a Poisson process (compare 
Eq. (0)). For Poisson distributed spheres with vanish- 
ing cumulants £,n(rn) = one recovers Af^ — 1. The 
specific Minkowski functionals depend on the radius r, 
the product density g, and all of the correlation func- 
tions ^„(xi,... ,x„). In particular, one obtains for the 
normalized Minkowski functionals in three dimensions 

iv = ^r^g) 



*l(r•;^?,Un}) = Mi e-'jA^o 



77M1 e-"^" 



$3(r; g, {U}) = - 377M1M2 + ^r,2M^ j e-"^" , 

(16) 

which can be compared to the analogous result (^ for 
Poisson distributed grains. The dependence on higher- 
order correlation functions enters only into a finite number 
of relevant coefficients M^,. The universal polynomial form 
of the mean values ( p^ is related to the additivity (^ of 
the Minkowski functionals due to the decomposition into 
specific terms M^,. We also define the deviations 



5M,{r) = (1-M,(r))^ 

r] 



(17) 
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from the uncorrelated Poisson values Af ^ (r) = 1 with 



0. 



We tried to use the specific Minkowski functionals M,^ 
or SMi, to compare the cluster distribution with our mod- 
els. However, due to the nonlinear dependence of the 
on the measured <i>,y, the relative errors are significantly 
enlarged, compared to the errors of the The discrimi- 
natory power of the $i^'s is lost. This may be understood 



^„ = for n > 2, defines the Gauss-Poisson point process: 

\ogG[h + 1] = g dxi /i(xi) + 

Jr'' 

+ ^ [ dxi / dx2 6(|X1-X2|) /l(xi)Mx2). (19) 



in det a|il by solving Eg. ( 1£) for the Mi^ and inserting an 
error $i,(r) — $iy(r) + A<I>^(r) for the measured values 
of the Minkowski functionals $1/. Expanding in powers of 
A$y one gets AAI^, (x e~^^'^" for any M^,. The errors 
of the specific functionals increase exponentially for large 
77. It is clear that the are the preferable choice in the 
comparison of data with the models, but for the analytical 
calculations the and 6M^ axe more appropriate. 

3.2. The Gauss-Poisson process 

The notion of a Gaussian random field is well understood 
in cosmology (e.g. Bardeen et al. 1986| ): considering the 
density contrast (5(x) = p[x.)/pH — 1 the two-point cor- 
relation function ^2(^) = ]E[5(0)(5(r)] together with the 
mean mass density pn specifies the statistical properties 
of the mass density field completely (E is the average over 
several realizations of the random field). The higher cor- 
relation functions all equal zero. In the following we 
will show how to construct a "Gaussian" point distribu- 
tion. A detailed discussion, examples, and extensions t o 
higher-order processes is presented in Kerscher (2001 ). 
The defining property of this Gauss-Poisson point process, 
similar to the Gaussian random field, is that the higher- 
order correlation functions of the point set vanish: ^„ — 
for n > 2. Due to the discrete nature, and the demand 
for a positive number density g, some constraints on the 
two-point correlation function ^2 as well as the number 
density g emerge. 

In general, a point process may be specified by its 
probability generating function al (p.g.fl.) G\h] where h 
are suitable functions (see e.g. Daley fc Vere-Jones 1988 
Sect. 7.4; G[h] = TZ[h] as defined by Balian fc Schaeffcr 
1989| ). The p.g.fl. is the point process analogue of the prob 



The Gauss-Po isson point process ( Newman 197C , Milne fc 



Westcott 1972 ) is stochastically fully specified by its two- 



ability generating function of a discrete random variable 
(Kendall fc Stuart, 1977). The expansion of G[h] in terms 
of the (connected) correlation functions ^„ reads: 



iogG[;i + i] = ^ 



n=l 



/ dxi---/ dx„ ^„(xi, . . . ,x„) /i(xi) • • • /i(x„). (18) 



If we truncate this expansion after n — 1 we arrive at the 
p.g.fl. of the Poisson process. A truncation after n = 2, i.e. 



point correlation function ^2(^) and the number density 
g. However, g and ^2(^) may not be chosen arbitrarily. 
Milne fc Westcott (1972 ) showed that the Gauss-Poisson 
process is only well-defined if two constraints are satisfied. 
In Appendix ^ we give a detailed derivation. 



A simplified version of the constraint (A.4) is 



g / dy 6(|y|) < 1. 



(20) 



This tells us that sitting on a point of the process on av- 
erage at most one other point in excess of Poisson dis- 
tributed points is allowed. The constraint (A. 5) implies 



^2{r) > for any r, 



(21) 



hence only clustering point distributions may be modeled 
as a Gauss-Poisson process. 

Clearly the question arises, what is wrong with the 
simple picture that we start with a Gaussian random field 
and "Poisson sample" it to obtain the desired point dis- 
tribution. The answer is that a Gaussian random field is 
an approximate model for a mass density field only if 
the fluctuations are significantly smaller than the mean 
mass density. Otherwise negative mass densities (i.e. nega- 
tive "probabilities" for the Poisson sampling) would occur. 
Only in the limit of vanishing fluctuations a Poisson sam- 
pled Gaussian random field becomes a permissible model. 
However, in this limit we are left with a pure Poisson pro- 
cess. 

3.3. Simulating the Gauss-Poisson process 



As discussed by Daley fc Vere-Jones (1988) any Gauss- 
Poisson process equals a rather simple type of Poisson 



cluster processes (for details see Kerscher 2001). A Poisson 
cluster processes is a two-stage point process. First we dis- 
tribute parent points y (the supercluster centers) accord- 
ing to a Poisson process with number density gp and then 
we attach to each parent a second point process (the su- 
percluster). In this specific example the supercluster con- 
sists only of one or two points with probability qi (y) and 
92 (y), respectively. We have qiiy) + 92(7) = 1 and the 
first point is the supercluster center itself. The probabil- 
ity density /(|x — y|) determines the distribution of the 
distance |x — y| of the second point x to the superclus- 
ter center with /dx/(|x|) = 1. The p.g.fl. of this Poisson 
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cluster process is given by 

logG[h+l] = / dxi gpil + q2 



(xi)) h{xi 



dxi / dx2 £ipq2(xi)/(|xi - X2I) /i(xi)/i(x2), 

(22) 

which equals the p.g.fl. (|l^) for the Gauss-Poi sson pro- 
cess fo r g = gp(l + 52) and 6(0 = 2^ /(r) ( paley fc 
Vere- J ones 1988 Sect. 8.3, Kerscher 200 1| ). Hence, every 



Gauss-Poisson process is a Poisson cluster process of the 
above type, and vice versa. For a Gauss-Poisson process 
it is necessary that the parent points, in our case the su- 
percluster centers, are distributed according to a Poisson 
process. Any deviation from a pure Poisson process, either 
in the direction of clustering or regularity, implies the pres- 
ence of higher-order correlation functions in the distribu- 
tion of the galaxy clusters. There are indications that the 
distribution of galaxies or galaxy clusters shows some reg- 



ularity on large scales (Broadhurst ct al., 1990; Einasto 
et al., 1997| ). Hence, an unambiguous detection of such 
large scale regularity in the upcoming large redshift sur- 
veys would also strengthen our findings of non-Gaussian 
features on large-scale (see below). However, regular fea- 
tures on large scales are not necessary for higher-order 
correlations to be present. 

We may generate realizations of the Gauss-Poisson 
process for a given number density g and two-point cor- 
relation function ^(r), fulfilling the constraints ( A. 4 ) and 
(P). With C2 = 4.dx 6(|x|) and /„,dx /(|x|) = 1 
we can calculate the quantities needed in the simula- 
tion: /(r) = 6(?')/C2, 92 = 2^=ffc' 91 = 1 - '?2, and 
9p = £*(! — ^^2/2). The constraint (|2(]| ) impHes C2g < 1. 
The simulation is carried out in two steps: 

— First we generate the parents (the supercluster cen- 
ters) according to a Poisson distribution with number 
density gp. 

— For each supercluster center y we draw a uniform ran- 
dom number (7 in [0, 1]. If (7 < qi, then we keep only 
the point y. If q > qi then additionally to the point 
y we chose a random direction on the unit sphere and 
a distance d with the probability density / and place 
the second point according to them. 

To get the correct point pattern inside the sample window, 
one also has to use supercluster centers outside the window 
to make sure that any possible secondary point is included 
in the sample inside the window. 

As an illustration of this procedure we calculate the 
two-point correlation function ^ from the sample L20. 
The 6 together with the number density satisfy the con- 
straint ( pO| ) (sec Table ^): gC2 < 1. We use this empirical 
6 as an input to the simulation algorithm outlined above. 
Fig. 1^ illustrates that these simulated Gauss-Poisson sets 
are indeed able to reproduce the observed two-point cor- 
relation function. Even the dip of ^ at 20/i~^Mpc, is well 
reproduced by the simulated point sets. By construction 




10 

r [Mpc/h] 

Fig. 7. The two-point correlation function ^(r) esti- 
mated from the L20 sample (solid line) is shown together 
with 6(r) estimated from 200 realizations of a Gauss- 
Poisson process with the same ^ and number density as 
for the observed clusters (short dashed line, shaded one-tr 
area). The long dashed hue is for ^(r) = (30/r)^-®'^. 

no higher-order correlations are present in the simulated 
point sets in the mean. 

3.4. Minkowski functionals of a Gauss-Poisson process 

In the following the Minkowski functionals for a Gauss- 
Poisson process will be given. Truncating after the second 
term in Eq. (|l^) one obtains the correlated average of the 
specific Minkowski functionals 



M.(r;g,6) =1 - 2''-'rj SMf^\r:,^2) 



(23) 



<5Mf^'(r;6) 



1 



dx6(|x|) /.(r;|x|) 



B2,.(0) 



with the volume Mo(i?r) of a sphere of radius r. From 
Eq. d23) and using the expressions for Iy{r\s) given in 
Eq. ( |ll| ) we can calculate the normalized Minkowski func- 
tionals of the Gauss-Poisson process according to Eq. (|l6|) . 
For example we obtain 



with 

Mo = 1 



^o(r-;g,6) = (1-e 
7] 47r 



-77M0 



2Mo{Br 



2r 



3 s 
22^ 



-(-)' 



and similar expressions for $1, $2, and $3. 
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The expansion to linear order in r] oc r^g allows 
us to describe the MFs only for small radii or low number 
densities. For a Gauss-Poisson process both and r/ ap- 
pear non-linear ly in Eg . (p^ ) via Eqs. |2^, |2^). Contrary 
to the approximation (|10D which is only valid for 77 <C 1, 
the Minkowski functional of the Gauss-Poisson process, 
are valid for all rj. A Gauss-Poisson process does not imply 
the linearity of the MFs in and 77. 

For a Poisson distribution with ^2(?') = one obtains 



i.e., one recovers Eq. (||). 
form 



, SAC\r)^0, (24) 
Assuming an algebraic scaling 



(25) 



for the centered two-point correlation function one ob- 
tains with 7 < 3 the general result 



A 



with the amplitudes 



^(7) 
^1(7) 

^2(7) 
^3(7) 



24 

3-7 
12 

3-7 

6 

3-7 
3 

3-7- 



36 
4 — 7 

12 
4-7 ' 

6 

4-7 



12 

6-7' 

Tl/2 



3 J dy y cos^ y sin 



2-7, 



(27) 



(with g 



In the limiting case 7 

the result ^0 = ^1 = ^3 = 1 and A2 = ^ 
averaged Minkowski functionals 



0) one recovers 
1 



^ for the 



dx 



M^{Br)MoiBr) 



(28) 



of the intersection of two overlapping spheres with dis- 
tance |x| and radius r (see kinematic formula (^)). 

One has to be careful in the interpretation of Eq. ( |27|) , 
since a scale-invariant two-point correlation function ( P5[ ) 
with 7 < 3 does not satisfy the constraint (|2^), as required 
for the existence of a Gauss-Poisson process. A cut-off, i.e. 
^2(5) — for s > Tc, has to be imposed below some ra- 
dius Vc to guarantee the constraint (pO|). The maximal al- 
lowed Tc is depending on the number density through (pO|). 
As can be seen directly from Eq. (23), the scaling behav- 
ior ( p6| ) as well as the amplitudes in Eq. (27) are still 
correct for 2r < Vc- Only on larger scales additional terms 
depending on the cut-off will emerge. In such a scale in- 
variant Gauss-Poisson process the specific MFs 6M^{r) 
should show the general scaling form ( p6[ ) which may be 
used to test for an algebraic two-point correlation function 

The actual measured specific MFs still depend on 
the density g and all of the correlation functions 
^("'^(xi,... ,x„). The deviations of the measured MFs 
from the expressions (^^ for a Gauss-Poisson process will 
be used as a measure for the relevance of higher-order cor- 
relations among the points (see Sect. 0). 



To facilitate future applications we will also quote the 
Minkowski functionals of a Gauss-Poisson process in two 
dimensions (Mecke, 1994). With rj ~ nr'^g the reduced 



Minkowski functionals (compare Eq. (16) for three dimen- 
sions) 



^2(r;e,{Cn}) 



-vMo 



Mie 

{M2-rjM\)e-'^'^\ 



(29) 



with the specific Minkowski functionals given by Eq. (|15|) . 
For a Gauss-Poisson process the specific Minkowski func- 
tionals are given by the Eq. (p3|), the volume Mo(-Br) = 
Trr^ of a disc of radius r, and the functions 



h[r;s) 
h{T\s) 



1 



arcsm 



(30) 



e(2r-s). 



2r 



(26) If the correlation function decays algebraicly ^(r) = 
(•so/'")^i ^ = |xi — X2I, with an scaling exponent 7 and 
a correlation length sq one finds for the intensities of the 
Minkowski measures ( p^ for homogeneously distributed 
discs Bj. of radius r in two dimensions the amplitudes 



(Mecke, 1994) 



^0(7) = 2^ 
A (7) 

^2(7)= 2^, 



4 

2-7 



2g(7) 
5(7), 



/(7), 



with the functions 



/(7) 



4 r(i-i) 
Ar(3-i) 



and 



ail) - J 




' arcsin(?/). 



(31) 



(32) 



(33) 



One also recovers the result = 1 in the limit 7-^0 
(with g ^ 0) for the averaged Minkowski functionals of 
the intersection of two discs defined in Eq. (|2|). In gen- 
eral a cut-off in the scale invariant correlation function is 
needed to make this model well-defined. See the comments 
above. 

4. Non-Gaussian morphology of the galaxy 
cluster distribution 

In this section we compare the Minkowski functionals de- 
termined from the cluster distribution with the Minkowski 



functionals of a Gauss-Poisson process. In Sect. 3.2 



showed that the number density g and the two-point cor- 
relation function ^2 ('') have to fulfill constraints in order to 
allow them to serve as the ingre dien ts for a Gauss-Poisson 
point process. The constraint (A^) implies £,2{r) > for 
all r. There are indications from the analysis of the flux- 
limited Reflex catalogue that the two-point correlation 
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function of the cluster distribution becomes negative on 
scales of 40-50/1" ^Mpc ( poUins et al., 2000] ). The viola- 
tion of constraint (A. 5) already tells us that the cluster 
distribution exhibits non-Gaussian features even on such 
large scales. Due to the limited number of clusters and 
the smaller extent of samples we can not detect this zero 
crossing unambiguously in the volume-limited samples we 
analyzed. To obtain a well-defined model of the two-point 
correlation function we impose a cut at rc = 48/i~^Mpc 
with (,2(r) = at r > Tc. 



Another more stringent constraint is Eq. (A.4) which 
may be cast into the form 



9 



dy 6(|y|) = gC2<i. 



(34) 



As can be seen from Table the smaller volume-limited 
samples L05-L12, with their higher number density clearly 
violate this constraint. Hence, already by inspecting the 
two-point correlation function ^2 in conjunction with the 
number density g we can conclude that there are non- 
negligible higher-order correlations in the point distribu- 
tion of galaxy cluster in the samples L05 and L12. 



In Sect 3.3 we showed how to simulate a Gauss-Poisson 



process for a given two-point correlation function ^2 and 
g. The violation of the constraint (^) prohibits the sim- 
ulation of a Gauss-Poisson process corresponding to the 
samples L05 and L12. However for the samples L20 and 
L30 the constraint (|3^) is satisfied (see Table ||) and we 
may generate realizations of a Gauss-Poisson process with 
the same number density and the same two-point cor- 
relation function, as estimated from these samples, (see 
Sect. 3^ and especially Fig. |^). Both constraints ( po|) and 
( ^ ) are only necessary conditions for the existence of 
a Gauss-Poisson process. Still higher-order correlations 
may be present in the cluster distribution. We calculate 
the Minkowski functionals of these Gauss-Poisson sam- 
ples and compare them with the Minkowski functionals 
of the observed cluster distribution. Since correlations of 
any order enter the Minkowski functionals (see Eq. (||) or 
Eqs. dis] ) and (16)), deviations of the Minkowski function- 
als of the cluster distribution from the Minkowski func- 
tionals of the Gauss-Poisson process indicate the pres- 
ence of higher-order correlations even in these deep sam- 
ples. To facilitate the comparison we use the normalized 
Minkowski functionals 



M^Br) g{ry 



(35) 



To get unbiased estimates of the MPs we employ a bound- 
ary corrected estimator where we shrink the observational 
window W{r) with increasing radius r of the spheres 
( iMecke et al 1994 |Schmalzing fc Diaferio 20001 ). As 
an estimate of the number density inside W{r) we use 
g{r) = N{r)/\W{r)\ with N{r) the number of points in- 
side the shrunken window with the volume |W^(r)|. 

In Fig. [s] the results of our comparison are shown, 
where the ^'"^ plotted against -q — £i(47r/3)r^. We 



used the empirical two-point correlation function to gen- 
erate the realizations of the Gauss-Poisson process inside 
the sample geometry of the Reflex cluster catalog. The 
shaded one-cr area with the short dashed line in the center 
was estimated from one hundred realizations. The initial 
slope of the MPs of the cluster distribution (solid line) is 
well approximated by the expression ( p^ . But already for 
fairly small 77 this linear approximation breaks down. As 
discussed in Sect. 3.4 a Gauss-Poisson process does not 



imply the linearity of the MPs in 77, which is readily ob- 
served in Pig. H It is necessary to compare the measured 
values of M^{r) with Eq. (^^ which is available in the 
analytic form only because of additivity. So our heuristic 
argument in the beginning (to look for additive integral 
information on higher correlations) turns out to be use- 
ful in deriving analytic results which are necessary for the 
comparison with measured values. Over the whole range of 
scales probed, the normalized volume $0 and surface area 
$1 of the cluster samples are consistent with the Gauss- 
Poisson process. However, both the normalized integral 
mean curvature $2 and the normalized Euler character- 
istic $3 are lowered with respect to the Gauss-Poisson 
process, clearly outside the one-cr range. The deviations 
are especially prominent for radii larger than 30/i~^Mpc. 
This is a firm indication that higher-order correlations 
are necessary to account for the shape and topology of 
the cluster distribution for such large scales, given by the 
radii of the spheres. 

In Fig. ^ also the MPs of a Gauss-Poisson process with 
a scale invariant two-point correlation function ^2(?') = 
(so/r)^ are shown (see Eqs. (M), (p7|) and (p^)). The expo- 



),(!!) 

bylC£ 



nent is 7 = 1.83 as determined by |Collins et al. (2000 ) and 
So = 30/i^^Mpc. These parameters give a reasonable fit 
to the wiggly two-point correlation function determined 
from volume-limited sample L20 (Fig. |^). The normalized 
volume $0i surface area $1 and integral mean curvature 
$2 of a Gauss-Poisson process with this scale invariant 
^2 follow closely the corresponding quantities determined 
from the Gauss-Poisson process using the two-point cor- 
relation function from the data. A significant difference 
between the two Gaussian models shows up only in the 
Euler characteristic $3. 

To quantify the deviation of the distribution of galaxy 
clusters from a Gauss-Poisson process we perform a 



non-paramet ric significance test ( Besag fc Diggle, 1977 ; 
^toyan, 2000| ). Using M — 10 equidistant radii in the 
range from 30/i~^Mpc to 45/i~^Mpc we define the "dis- 
tance" 



1 



(36) 



between the normalized Minkowski functionals <I>J; of the 
sample labeled with k and the mean value <I>^^ of the 
Gauss-Poisson process calculated from one hundred real- 
izations inside the same sample geometry as the Reflex 
samples and using the empirical g and ^2- Additional to 
the distances of the cluster sample L20, we also cal- 
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Fig. 8. The reduced Minkowski functionals <i>^ of the volume-Hmited sample L20 (soUd hnc, see also Fig. ||) compared 
to the Minkowski functionals of a Gauss-Poisson process (short dashed line, shaded one-a area) with the empirical 
two-point correlation function as input (see Fig. ^). The long dashed line is obtained from the scale invariant model 
(see Eq. (|2|)) with sq = SO/i^^Mpc, 7 = 1.83 and the amplitudes Aq = 6.8, Ai = 4.73, A2 = 2.70, and A3 = 2.56 
according to Eq. (|2^). The mean Minkowski functionals of a Poisson process are shown as a dotted line. The errors 
for the Poisson process (not shown) are smaller than the errors from the Gauss-Poisson process. 



culate the distances d'^ for k — 1 ... 99 realizations of 
the Gauss-Poisson process (these 99 realizations are inde- 
pendent from the realizations used to calculate the mean 
$^^). We order these distances including d^^'^ ascending. 
If dl;^'-' is under the five highest distance values, we may ex- 
clude a Gauss-Poisso n process with a significance of 95% 
(see the comments by Marriott 1978 concerning the signif- 
icance level). The beauty of this Monte-Carlo significance 
test is that we neither make assumptions about the distri- 
bution of the Reflex clusters, nor about the distribution 
of the errors of the MFs in the model. 

In Table ^ the rank of the cluster sample L20 within 
the ordered list of distances is given. As expected from 



Table 3. The rank of the cluster sample L20 within the 
ascending list of one hundred ordered distances dj; . A rank 
larger than 95 indicates rejection with significance of 95%. 



V 





1 


2 


3 


rank 


26 


52 


97 


96 



the visual impression in Fig. g the volume $0 ^'^^ the 
surface area $1 are consistent with the expectation from 
a Gauss-Poisson process. Both the integral mean curva- 
ture <I>2 and the Euler characteristic $3 allows us to reject 
the hypotheses that the cluster distribution stems from a 
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Gauss-Poisson process at a significance level of 95%. This 
result is stable against extending or shrinking the radial 
range. One may also use the overall mean number density 
Q of the full sample in Eq. ( |35| ) instead of e(r). 

We implemented this Monte-Carlo test using the em- 
pirical two-point correlation function as input to the simu- 
lations of the Gauss-Poisson process. For a scale-invariant 
C2('") = (30/i^^Mpc/r)^-*'^ the normalized Euler charac- 
teristic $3 according to Eq. ( ]26| ) seems to be in agree- 
ment with the observed MPs. However, the family of MPs 
provides us with a consistency check, still the integral 
mean curvature $2 from the scale invariant model and 
the data are differing, illustrating the relevance of higher- 
order correlations. Hence, also a Gauss-Poisson process 
with a scale-invariant correlation function is inconsistent 
with the data. 

We conducted a similar analysis for the cluster sample 
L30 (see Pig. H). Again $0 and <i>i fall within the one-cr 
range of the Gauss-Poisson process. The $2 an d $3 only 
marginally stand out. As discussed in Sect. 2.3 this sam- 
ple hardly allows for a discrimination from the Poisson 
process, which can be explained by Eq. ( [l0| ) and the sig- 
nificantly lowered number density. Nevertheless the same 
tendency can be observed as for the sample L20 although 
the statistics does not allow a discrimination. 

5. Summary 

The Reflex cluster catalogue is well suited for study- 
ing the large-scale structure of the Universe. The detec- 
tion of the clusters is based on their X-ray flux, allow- 



are needed to allow for the increased variance. Clearly, 
the relevance of these higher-order correlations has to be 
checked independently, e.g. using the MPs. Due to the 
decreasing number density of galaxy clusters the deeper 
volume-limited samples L20 and L30 comply with the 
constraints. A Gauss-Poisson process based on the ob- 
served correlation function becomes feasible as a model. 
MPs summarize the influence of the two-point correlations 
and higher-order correlations on the morphology of large- 
scale structure. They include correlations of any order in 
an integral way. We calculated the MPs for a general corre- 
lated point set. Detailed results were given for the Gauss- 
Poisson process. To quantify higher-order correlations in 
the cluster distribution we compare the analytical known 
MPs known for the Gauss-Poisson process with the ac- 
tual observed MPs of the cluster distribution. Two of the 
four MPs, the volume and the surface area, are consis- 
tent with the Gaussian model. However a clear detection 
of non-Gaussian features at large scales was possible with 
the integral mean curvature and the Euler characteristic. 

The definition of the Gauss-Poisson process directly 
lead to a method for simulating Gaussian point distri- 
butions. With such simulated point distributions we per- 
formed a non-parametric Monte-Carlo test. The main re- 
sult is that we can exclude a Gauss-Poisson process as a 
viable model for the distribution of galaxy clusters at the 
significance level of 95%. 

Non-Gaussian features seen in the distribution of 
galaxy clusters may be already imprinted on the initial 



ing the construction of a flux -limited sample. X-ray se- 



density field (see e.g. Linde & Mukhanov 1997), or may 
be a result of topological defects (see e.g. shellard & 
Brandenberger 198S). We would like to point out that 



lected cluster catalogues are not impaired by projection 
effects. Moreover, the flux-limit, together with the well 
documented selection effects allows the extraction of clean 
volume-limited samples. 

We calculated Minkowski functionals of a series of 
volume-limited samples, extracted from the Reflex clus- 
ter catalogue. The comparison with the MPs of Poisson 
distributed points revealed similar features as detected 
in the Abell/ACO cluster sample ( Kerscher et al., 1997 ). 
Although the number of clusters in the samples is always 
less than one hundred, MPs allow for a sensitive and dis- 
criminatory analysis. The stability of the results obtained 
from this small number of points can be attributed to the 
additivity property of the MPs, which served as a con- 
struction principle. 

Our aim was the quantification of non-Gaussian fea- 
tures in the large-scale distribution of clusters, therefore 
we first gave a precise definition of a Gaussian point 
distribution, the Gauss-Poisson process. Contrary to a 
Gaussian random field, constraints for number density and 
the two-point correlation function arise. In the smaller 
volume-limited samples L05 and L12 these constraints 
are violated. Hence, a Gauss-Poisson process with the 
observed density and two-point correlation function does 
not exist. This is an indirect detection of higher-order 
correlation functions. Higher-order correlation functions 



also explanations facilitating Gaussian initial conditions 
are possible. Introducing a threshol d and considering only 
peaks in a Gaussian density field Bardeen et al. (1986 ) 
could show that the point distribution of the peaks has 
non zero higher-order correlations ^„ 7^ for n > 2. Still 
the importance of the higher-order correlations on large 
scales comes as a surprise within this model. On physical 
grounds, the peak biasing picture may only serve as a first 
approximation. Evolving a Gaussian density field in time 
using the linear approximation leads to larger and larger 
regions with a non-physical negative mass density. Only 
the non-linear evolution of the density field can remedy 
these shortcomings, allowing on the one hand for high den- 
sity peaks with an over-density of several hundreds, and 
on the other hand allowing voids with a negative density 
contrast always larger than minus one. At the peaks of 
this non-linear evolved density fields one may assume the 
clusters to reside. As already discussed in the introduc- 
tion non-Gaussian features in the large-scale distribution 
of mass, like walls and filaments, are predicted both by 
the Zel'dovich and related approximations as well as by 
A^-body simulations, both based on Gaussian initial con- 
ditions. These structures in the mass distribution, perhaps 
amplified by a biasing mechanism, can be associated with 
the non-Gaussian structures observed in the large-scale 
distribution of Reflex galaxy clusters. Our results sug- 
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Fig. 9. The reduced Minkowski functionals of the volume-Hmited sample L30 (see also Fig. The same conven- 
tions as in Fig. H apply. 



gest that within these scenarios, using Gaussian initial 
conditions, it is necessary to consider non-linear models 
to describe the observed large-scale structures. 
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With Ai as an infinitesimal volume element centered 
on the origin and Aj equal to some volume A the first 
constraint (A.4) implies the simplified constraint (pO[). 
Considering two volume elements Ai = dV^ and Aj ~ dVj , 
then Eq. ( |AT5| ) impHcs Eq. @. 



Appendix A: The constraints 

In this appendix we recall the derivation of the constraints 
on the number density g and ^2(?") for the Gaus s^Poisson 
process as pres ented by Milne fc Westcott (1972| ) (see also 
Kerscher 20"0ll ). 

Consider k compact disjoint sets Aj, and let rij = 
N{Aj) be the number of points inside Aj. The probability 
generating function of the fc-dimensional random vector 
(ni, • • • , Hk) is then 



n 



(A.l) 



Together with a continuity requirement the knowledge of 
all finite-dimensional probability generating functions Pk 
determines the p.g.fl. G [h] (see Eq. ([T^ )) and the point 
process completely (e.g. [Westcott 1970| ). Setting 



Mx)- 1-^(1 -z,)1a,(x), (A.2) 
i=i 

one obtains the probability generating function of these 
finite-dimensional distributions from the probability gen- 
erating functional of the point process: Pfc(z) — G[h]. 
Here ]lyi(x) is the indicator-function of the set A, with 
flA(x) = 1 for X E A and zero otherwise. 

Since Pfc(z) is a probability generating function of a 
random vector, it is positive and monotonically increasing 
in each component Zi, hence logPfc(z) is non-decreasing. 
Inserting Eq. ( |A.2| ) into the probability generating func- 
tional of the Gauss-Poisson process (|l^) one immediately 
obtains 



dzi 



q\Ai\ + 



k 

+ g'J2f dx / dy6(|x-y|)(z, -l)>0 (A.3) 

j^l^Ai J Aj 

for any Zj > 0, where \Ai\ is the volume of the set Ai. We 
chose Zj = 1 for all j ^ i and set the remaining Zi either 
to = or ^ 1. Then the following two constraints 
emerge: 

^ ^ dx|^ dy6(|x-y|)<l, (A.4) 
/ dx / dy6(|x-y|) >0, (A.5) 

J A, J A, 



for any subset A,,Aj of M^. [Milne fc Westcott (1972|) 
showed that these two conditions are necessary and suffi- 
cient for the existence of the Gauss-Poisson process. 



